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METHOD ^F^FINE SYNC HRONIZATION T O— A SXGNALi_RE_CEIVED^FROM 

A TRANSMISSION CHANNEL 



The present invention relates to a method of fin 
synchronization to a receive signal corresponding to 
5 reference signal transmitted in a transmission channe 
In a transmission system, in particular a radio 
transmission system, a transmitter transmits a reference 
signal to a receiver in a transmission channel. One of 
the first operations that the receiver must perform is to 
10 synchronize to the receive signal. This problem is well- 
known to the skilled person and it is therefore 
unnecessary to remind the reader here of the various 
techniques that have been used to obtain such 
: - ; synchronization . 

*jj 15 If the signals are digital signals, it is usual to 

y evaluate the synchronization error by means of a time 

;Ts unit which is the time difference between two successive 

bits of a signal and which is referred to as the bit 
^ § period. It appears that the available prior art 

20 solutions cannot achieve synchronization to an accuracy 
"2 much better than one bit period. 

3 This accuracy can prove insufficient in some cases. 

This is because synchronization gives the transit time of 
the signal in the transmission channel, in other words 

25 the transmission time between the transmitter and the 
receiver. This data is important in a duplex 
radiocommunicat ions system in which a base station 
communicates with a terminal because the base station and 
the terminal have respective transceivers and the 

3 0 terminal must operate in such a way that the signal it 

transmits arrives at a precise time with reference to the 
clock of the base station. To achieve this, it is 
naturally necessary for the terminal to know the time 
taken for the signal to reach the base station. 

35 The transmission time is also a direct reflection of 

the distance between the transmitter and the receiver. 
Clearly, the accuracy of this distance is of fundamental 




importance when it is a question of locating the terminal 
by identifying its position relative to one or more base 
stations. The general problem of locating a terminal is 
a very current concern because of its applications, which 
include cell handover strategies in cellular networks, 
for example. The security field should also be 
mentioned, whether in connection with identifying the 
geographical location of the source of an emergency call 
or the position of a stolen vehicle equipped with the 
terminal . 

An object of the present invention is a 
synchronization method whose accuracy is much better than 
one bit period. 

In accordance with the invention, the method of fine 
synchronization to a receive signal corresponding to a 
reference signal transmitted in a transmission channel 
includes the following steps: 

- selecting a source signal producing a characterization 
signal after it has passed through said transmission 
channel , 

- establishing a characterization matrix for estimating 
the covariance of the characterization signal, 

- identifying dominant eigenvalues which are the highest 
eigenvalues of the characterization matrix, 

- calculating the correlation function of the source 
signal with the sum of the eigenvectors associated with 
the dominant eigenvalues, and 

- searching for the first maximum of the correlation 
function. 

If the time increment adopted for calculating the 
correlation function is made sufficiently small, and in 
all cases much less than one bit period, the above method 
achieves very good accuracy. 

A first option is for the number of dominant 
eigenvalues to be predetermined. This number typically 
represents 20% to 30% of the dimension of the 
characterization matrix . 
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A second option is for the ratio of the sum of the 
dominant eigenvalues to the sum of the all the 
eigenvalues to be greater than or equal to a 
predetermined percentage. In this case the percentage 
5 adopted is often greater than 90%, for example 95%. 

A third option is for the method also to include a 
step of estimating the additive noise in the transmission 
channel, the dominant eigenvalues being such that their 
sum is less than or equal to the sum of all the 
10 eigenvalues less the additive noise. 

Also, the additive noise is estimated by normalizing 
the instantaneous noise, which is evaluated from the 
*Jf receive signal, the reference signal and an estimate of 

oi the impulse response of the transmission channel . 

s "*f 15 The expression for the instantaneous noise is 

Q advantageously: N 0 = S - A.X, where A is the transmission 

matrix associated with the reference signal . 

Whichever of the above options is adopted, the 
lj characterization matrix results from a smoothing 

I J 2 0 operation. 

M= In a preferred embodiment, the characterization 

;rf signal is an estimate of the impulse response of the 

transmission channel . 

The characterization signal can instead be the 
25 receive signal. 

The present invention emerges in more detail from 
the following description of proposed embodiments of the 
invention, which is given by way of illustrative example 
and with reference to the accompanying figures, in which: 
30 - Figure 1 shows a first variant of the invention, and 
- Figure 2 shows a second variant . 

Elements common to both figures are allocated the 
same reference numbers . 

The receiver has already achieved coarse 
35 synchronization to the receive signal, to an accuracy of 
the order of one bit period, using any of the available 
prior art solutions. 



The receive signal corresponds to a reference signal 
produced by the transmitter and known to the receiver. 
The reference signal can be known a priori, i.e. it can 
be a training sequence made up of identified symbols. It 
can also be determined a posteriori using techniques 
referred to generically as blind probing. In this case, 
during the synchronization procedure, the receiver 
regenerates from the receive signal the series of symbols 
forming the reference signal. 

It is first necessary to characterize transmission 
between the transmitter and the receiver. To this end, a 
source signal produced by the transmitter is selected 
which yields a characterization signal in the receiver 
after transmission over the channel . 

Of course, if the source signal is the reference 
signal, the characterization signal is the receive signal 
itself. This is not always the optimum solution, 
however, in terms of the complexity and performance of 
the method of the invention. 

Another solution is to retain a modulated pulse as 
the source signal, in which case the characterization 
signal is the impulse response of the transmission 
channel . 

For example, the GSM digital cellular mobile 
telephone system uses a training sequence made up of 2 6 
symbols, the impulse response being generally estimated 
with five coefficients since it is accepted that the 
dispersion of the channel is equal to 4 . 

In this case, the receive signal has a maximum 
dimension of 22, which is significantly larger than that 
of the impulse response. 

Two embodiments of the invention are therefore 
examined in succession and with reference to Figure 1, 
starting with the situation in which the source signal is 
a modulated pulse. 

Estimating the impulse response does not give rise 
to any problems in itself because many methods are 
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available for this, for example the least squares method, 
which is described in particular in patent applications 
FR 2 696 604 and EP 0 564 849. Briefly, that technique 
uses a measurement matrix A constructed from the training 
sequence TS of length n. The matrix A has (n-d) rows and 
(d+1) columns, where d represents the dispersion of the 
channel. The element in the ith row and the jth column 
is the (d+i-j)th symbol of the training sequence. a.± 
denotes the ith symbol of a sequence TS of 26 symbols: 

a 4 a 3 a 2 a x a 0 
a 5 a 4 a 3 a 2 &i 

^■6 a s a 4 a 3 a 2 
A = a 7 



a 25 a 21 

The training sequence is chosen so that the matrix 
A t A can be inverted, where the operator . fc represents 

2 0 transposition . 

In the receive signal S, the first four symbols s 0 
to s 3 are ignored because, given that the dispersion of 
the channel is 4, they also depend on the unknown symbol 
transmitted before the training sequence. The receive 
25 signal is defined hereinafter as a vector S having for 
its components the received symbols s 4 , s 5 , s 6 , s 25 . 

The estimated impulse response X therefore takes the 

form: 

X = (A fc A) _1 A fc .S 

3 0 The next step of the method of the invention is to 

establish a statistic of this impulse response, where 
"statistic" means a data set reflecting the average value 
of the impulse response over an analysis period. 

A matrix is therefore constructed for smoothing the 
35 various estimates X obtained during the analysis period 
to obtain an estimate of the covariance associated with 
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the impulse response. In this context the term 
"smoothing" is to be understood very broadly, i.e. as 
referring to any operation for smoothing or averaging the 
impulse response over the analysis period. 



A first example of smoothing entails calculating the 
average of the matrix XX h over the analysis period, which 
is assumed to include m training sequences: 

1 m 

L(XX h ) = — ^XX h 



m 



i 



The operator , h represents the Hermit ian 
transformation or complex conjugate transpose. 

A second example of smoothing entails, after the ith 
training sequence is received, updating the smoothing 
15 matrix (XX h ) obtained on the (i-l)th training sequence 

by means of multiplier coefficient a, this factor 
generally being known as the smoothing forget factor and 
ranging from 0 to 1 : 

Li(XX h ) = ocXiXi 11 + (1-ot) Li.! (XX h ) 
20 Any means can be used for this initialization, 

including the first estimate X obtained, or an average 
obtained as described for a small number of training 
sequences . 

For simplicity, the smoothing matrix L (XX h ) , which 
25 is in fact a statistical characterization matrix, is 
denoted L below. 

The method then includes a step of looking for 
(eigenvalue, eigenvector) pairs of the characterization 
matrix . 

3 0 This step is not described in more detail because it 

is well-known to the skilled person. 

The eigenvalues k ± are then classified in decreasing 
order. Their sum corresponds to the energy of the 
characterization signal X made up in part of a wanted 
35 signal which images the source signal and in part of the 
additive noise N of the transmission channel. 

The dominant eigenvalues, the ones which are the 
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highest, represent the wanted signal, whereas the lowest 
eigenvalues represent noise. 

A first option is to retain a predetermined number 
of dominant eigenvalues. For example, the first two 
5 eigenvalues X x and X 2 are retained for an impulse response 
with five coefficients. 

A second option is to consider the wanted signal to 
have an energy which is a predetermined fraction f. of the 
energy of the characterization signal. Thus using the 
10 notation X ± for the eigenvalues with i. varying from 1 to 
P, there will be d dominant eigenvalues, where d is 
obtained as follows: 

d d+l 

^— <f and >f 
p p 

i=l i=l 

15 

The fraction f can be fixed a priori, for example at 
95%. This fraction can also be derived from the signal- 
to-noise ratio of the receive signal obtained elsewhere. 

A third option, and doubtless that offering the best 
2 0 performance, is for the additive noise N to be estimated 
directly from the receive signal and from the measurement 
matrix A. If N 0 denotes the noise vector affecting the 
receive signal : 

S = AX + N 0 

2 5 Given that the vectors S and N 0 have 2 2 components, 

the additive noise N can be expressed as follows: 

N = (^)(S-AX) h (S-AX) 

3 0 This estimate of the additive noise can naturally be 

averaged and smoothed. 

Using the previous notation, and standardizing the 
energies : 
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x ; >N and j\j < N 
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Thus the dominant eigenvalues are obtained from a 
direct estimate of the noise. 

Whatever option is previously adopted, the next step 
of the method consists of calculating the correlation 
function of the source signal with the sum of the 
eigenvectors v ± associated with the dominant eigenvalues 



The source signal is oversampled relative to the bit 
period and is therefore denoted g(t) where t. represents 
time and is a discrete variable whose quantization 
increment is 1/32 bit period, for example. It is 
represented by a vector with the same number of 
dimensions as the characterization signal, i.e. five 
dimensions in this example. The correction function c(t) 
is calculated for t varying from -1 to +1 bit period, for 
example, using the following expression: 



The period between the source signal g(t) and the 
eigenvector v ± represents a scalar product, in the 
conventional way. 

The final step of the method looks for the value t 0 
of t that is closest to zero, which corresponds to the 
first relative maximum of the correlation c(t) . It is 
this particular value t 0 which gives the required 
synchronization error relative to the reception signal. 

The following complementary function c' (t) can also 
be considered: 



Note that the particular value t 0 mentioned above 
can also be obtained by seeking the value of t. closest to 




i = d + l 



zero, which corresponds to the first relative minimum of 
the complementary formation c* (t) . 

These two methods of obtaining the synchronization 
t 0 are therefore equivalent . 

Turning to Figure 2, a second embodiment of the 
invention is considered in which the characterization 
signal is the receive signal S, so that the source signal 
is now the reference signal, or, in the case of the GSM, 
the training sequence TS when it is GMSK (Gaussian 
minimum shift keying) modulated. 

The statistic of the characterization signal is 
therefore estimated by means of a characterization matrix 
which is now obtained by smoothing the various 
occurrences of the receive signal S. Once again, the 
term "smoothing" is to be understood in a very broad 
sense in this context. 

The characterization matrix L therefore takes the 
following form: 



The p' eigenvalues X* ± of the matrix are then looked 
for, and as in the first embodiment the d' dominant 
eigenvalues are identified. 

The correlation function f (t) of the modulated 
training sequence and the sum of the eigenvectors v'i 
associated with the dominant eigenvalues X* ± are now 
calculated . 

Once again, the training sequence g'(t) is 
oversampled and is represented by a vector with 22 
components. The correlation function therefore becomes: 




or 



L ± (SS h ) = aSiSi h + (l-cOLi^CSS 11 ) 



f(t)=i>(t).v\ 



1 = 



As before, a new complementary function f ' (t) can be 
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defined : 



p' 



f(t)= £g'(t).Vi 
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The method terminates in the same manner by seeking 



the first maximum of the correlation function f (t) or the 
first minimum of the complementary function f ' (t) . 

The invention can therefore be implemented in 
various ways, the essential point being to have access to 
10 a source signal and of the result of transmitting it, 
i.e. the characterization signal. 



